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Abstract High precision Superconductivity Gravity Gradiometers (SGG) are 
powerful tools for relativistic experiments. In this paper, we work out the tidal 
signals in non-dynamical Chern-Simons modified gravity, which could be mea¬ 
sured by orbiting SGGs around Earth. We find that, with proper orientations 
of multi-axes SGGs, the tidal signals from the Chern-Simons modification can 
be isolated in the combined data of different axes. Furthermore, for three-axes 
SGGs, such combined data is the trace of the total tidal matrix, which is 
invariant under the rotations of SGG axes and thus free from axis pointing 
errors. Following nearly circular orbits, the tests of the parity-violating Chern- 
Simons modification and the measurements of the gravitomagnetic sector in 
parity-conserving metric theories can be carried out independently in the same 
time. A first step analysis on noise sources is also included. 


Keywords Chern-Simons Theory • Models of Quantum Gravity • Experi¬ 
mental Relativity • Gravity Gradiometer 

1 Introduction 

Among the modified gravitational theories, the extensions to the Einstein- 
Hilbert action with second order curvature terms are of particular interest, 
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which may arise in the full, but still lacking, quantum theory of gravity as high 
energy corrections to GR, see [T]. The string theory inspired Chern-Simons 
(CS) modified gravity PMlffllMS] , with the additions of a parity-violating term 
R *R and a coupling scalar field 9, is one of such extensions of GR. Being a 
promising model, CS modified gravity has found connections with different 
fields such as gravitational physics, particle physics, string theory, loop quan¬ 
tum gravity, and cosmology, please consult [6] for detailed discussions. 

CS modified gravity now contains two classes of formulations, the non- 
dynamical and dynamical formulations, which are in fact two distinct theories. 
In the non-dynamical formulation, the CS scalar 8 is externally prescribed. The 
so called canonical choice is to set 9 as a linear function of coordinate time 
proportional inversely to a mass scale Mcs IS- While, in the more realistic 
but complicate dynamical formulation, the evolution of the CS scalar is then 
sourced by the spacetime curvature. The non-dynamical formulation serves 
now as a useful model that provides us insights into parity-violating theories 
of gravity. Up to now, the experimental tests and constraints on CS gravity are 
all based on observations from astrophysics and space based experiments. The 
first but weak bound on the canonical CS scalar 9 was obtained in [7] based on 
the results from LAGEOS I, II [5][Dj and Gravity Probe-B [ID| missions, which 
had constrained the CS mass scale to Mcs > 2 x 10~ 13 eU. The strongest 
bounds on the canonical scalar up to now was based on the data from double 
binary pulsars HU, which had the constraint Mcs > 4.7 x 10 10 eU as been 
revised in [12 . For dynamical CS gravity, the vacuum solutions outside the 
rotating black holes and stars were studied in the slow rotation approximation 
in dMHHS], and their possible tests can be found in mnm- Moreover, 
the parity-violating term R *R also leaves distinguishable signatures in gravi¬ 
tational waves, which may be captured by ground based or future space borne 
gravitational wave antennas, please consult HSJEOEE] for details. 

In experimental relativity, Braginskii and Polnarev had obtained for rela¬ 
tivistic gravitational theories an interesting spin-quadrupole coupling between 
rotating sources and orbiting quadrupole oscillators mm- Along this line, 
the principles of detecting Earth gravitomagnetic field with orbiting gradiome- 
ters are studied in a series of works of Paik, Mashhoon and their collabora¬ 
tors Today, gradiometers have already been employed in 

space based experiment such as the one on GOCE satellite m- The high pre¬ 
cision Superconductive Gravity Gradiometer (SGG) have been developed at 
the University of Maryland by Moody, Paik and their colleagues, and their 
performance level has already reached 2 x 10" 11 s~ 2 /Hz^ in 1990s |.'i2| . Re¬ 
cently, with the improvement in cryogenics and the magnetically levitated test 
masses, a new design scheme with 3 orders of magnitude improvement in the 
SGG performance over the frequency band 0.5 mHz ~ 0.1 Hz has been devel¬ 
oped by this group mm- It is also pointed out in [35; that a performance 
level of 10~ 15 s~ 2 /Hz? over this band is within the capability of the SGGs 
under development. One can consult for the detailed error analysis 

of the relativistic gradient measurements with SGGs in space. On the other 
hand, as pointed out by Alexander, Yunes and etc. j37j|[38] that, within the 
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weak field and slow motion limits, the non-dynamical CS gravity differs from 
General Relativity (GR) only in the gravitomagnetic sector. Based on these 
results, we study in this work the theoretical principles of testing the parity- 
violating non-dynamical CS gravity with SGGs in space. 


2 Signatures of gradient measurements in the non-dynamical 
Chern-Simons gravity 

This work is heavily based on the mission concepts studied in [27 p g p 5 ] , that 
the gradient force from Earth gravitomagnetic field are measured by orbiting 
SGGs along nearly circular orbits. We worked out in this section the signatures 
of the gradient observable in non-dynamical CS gravity. 


2.1 The Non-dynamical Chern-Simons modified gravity 


We first give a brief introduction to the non-dynamical formulation of CS mod¬ 
ified gravity, for detailed discussions please consult (S1I571I5B] . The geometric 
units c = G = 1 are adopted in the followings. The action for non-dynamical 
CS gravity reads 

S ■— Sqr T Scs T 


where 


SoR= lk I dix ^9 R , Scs = j J d*xV^gOR*R } ( 1 ) 


and S m att is the action from the matter fields that is independent of 9. g is 
the determinant of the metric and the Pontryagin density 

R*R=\e cdefRa bef R b acd . (2) 

The CS coupling field 9 is externally prescribed and depends on the specific 
theory that under consideration. 9 can be viewed as the deformation function, 
and the difference between CS gravity and GR is proportional to the defor¬ 
mation parameters V a 0 and V a Vb$. In the so called canonical CS coupling 9 
is a spatially isotropic function and depends linearly on the coordinate time t 
[5], therefore the deformation parameter contains only 9. 

The field equation of the non-dynamical CS gravity is obtained by varying 
the action with respect to the metric 

Rab - T^gabR + 16naC ab = 8nT ab , (3) 


where C ab is the 4-dimensional generalization of the Cotton-York tensor 


C ab = \/ c 9e cde{a \/ e R b) d + ^\7 c \7 d 9e efd ( a R b)c f 


f e ' 


(4) 
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The introduction of the new scalar degree of freedom 0 also gives rise to the 
new constraint 

\7 a C ab = ~l\7 b d(*RR) = 0. (5) 

8 

If the above constraint is satisfied, from Eq. (|3]) the Bianchi identities and the 
equations of motion for matter fields V a T ab = 0 are recovered, which rank the 
non-dynamical CS gravity a metric theory. 

In the weak field and slow motion limits, the Parametrized Post-Newtonian 
(PPN) metric of the non-dynamical CS gravity outside a compact source was 
carefully worked out in [37l|38]. As mentioned before, the non-dynamical CS 
gravity differs from GR only in the gravitomagnetic sector 

Sof = + x(rVx V),, (6) 

here r denotes the distance to the mass center of the compact source and Vi is 
the PN potential, see Appendix A The dimensionless parameter \ = 32 ^ afl is 
the new PN parameter for non-dynamical CS gravity, and the CS mass scale 
is defined as |37|l38] 

1 4 

Mcs = -- = —. (7) 

8nad X r 


2.2 The basic settings 

We restrict ourselves to the so-called “semi-conservative” metric theories, which 
are based on action principles and respect the conservation law of 4-momentum 
[3]J]. Therefore, we work with the metric theories with four relevant standard 
PPN parameters {7, /?, £, 07, a 2 } together with the additional CS parameter 
X, please see [3jEj| or Appendix [A] for the PPN formalism. The PN coordinates 
system {t,x 1 } outside Earth is chosen as follows. The mass center of Earth 
is set at the origin. The basis (gfs-) 0 is set to parallel to the direction of the 
Earth angular momentum J, (gfr)° is pointing to a reference star T and (572)“ 
determined by the right-hand rule (gfr) a x (gfj) a = (gf?) 0 - Such coordinate 
directions are tied to the remote stars, and the time t is measured by the 
observers at asymptotically flat region. Within our coordinate system the PN 
metric outside Earth reads 

3oo = -1 + 2U - 2/3U 2 - 2 & w + (2 7 + 2 - 2f )#! 

+2(3 7 - 2/3 + 1 + + 203 + 2(3 7 - 2^)04 

+2£_4 — (ay — 0 L 2 )w 2 U — a 2 W l w^Uij — 2ay w l Vi + 0(e 6 ), 

goi = -^(47 + 3 + ai - a 2 - 2£)Vi - ^(1 + a 2 + 2^)Wi 

+xr(v x vy - i(ai - 2a 2 )w l U - a 2 w 3 U lJ + 0(e 5 ), 

9ij = (1 + 2jU)6ij + C>(e 4 ), 
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please see AppendixjAjfor the PN potentials. For low and medium Earth orbits 
experiments, the magnitude of e is about 10 -5 . 

We model Earth as an ideal and uniform rotating spherical body, and for 
the effects from Earth gravity multiples on the gravitomeganetic field measure¬ 
ments with orbiting SGGs, please consult [5B]. The preferred-frame and the 
preferred-location effects are tightly constrained by observations, and we now 
have the upper bounds of the related PN parameters as a\ ~ 4 x 10 -5 , 012 ~ 
2 x 10 —9 , 03 ~ 4 x 10 -20 and £ ~ 10 -9 , please see Tab[l]or [39] for more 
details. Generally, the coordinate velocity w of the PPN coordinate system 
relative to the mean rest-frame of the universe is generally believed to be 
small, that w ~ 0(e) |40U411155] , Therefore, the gradients produced by the 
preferred-frame and the preferred-location effects in orbiting SGGs will be 
smaller than 10 - 21 s -2 , which is too small to be seen by the present day SGGs 
and will be ignored in this work. The above metric can then be cast into a 
rather simple form 


j _|_ 2 M _ 20 M 2 


(^r + 


2 r 4 




2 r 4 


-)J~ 


t , AXX_ T 

V r 3 ' 2 r 4 

( Ax 1 , 3\x 2 x 3 \ r 
V r 3 ' 2 r 4 


\ X\{x 1 ) 2 + (x 2 ) 2 - 2 (x 3 ) 2 \ j 
\ 2 r 4 J 


1 i grM 

' r 

o 

o 


0 

i i g iM. 

I r 

0 


xKx'y+i^V- 2 (x 3 ) 2 } 

2 r 4 

0 


0 

l _j_ 27 M 


J\ 


27 M i 

r> / 


( 8 ) 


where r = sjSijx 1 xi , A = 1 + 7 + |ai, and M, J are the asymptotically 
measured total mass and angular momentum of Earth 


M = 

J = 


[ p[ 1 + (7 + l)v 2 + (3y -2/3 + 1 )U + — + 3~/-}d 3 x, 
J P P 

J jo(x x v)d 3 x. 


For a satellite orbiting Earth with velocity v. one has the basic order relations 






(9) 


2.3 Gravity gradient signals in non-dynamical Chern-Simons Gravity 

For clearness, we set the spacecraft (S/C) that carrying the multi-axes SGG 
to follow a circular orbit 

a; 1 =acos>F, x 2 = a cos i sin !F, a ; 3 = a sin * sin <F, ( 10 ) 
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Fig. 1 The orbits under consideration are circular orbits with the longitude of ascending 
node ft = 0. As illustrated in the figure, the local tetrad carried by the S/C is defined as 
follows, (ei) a is set along the direction of motion of the S/C, (& 2 ) a along the radial direction 
r and (£ 3 )“ = (ei) a x (e 2 ) a . 


where the longitude of ascending node 17, the eccentricity e are set to be zero, 
i and a denote the inclination and semi-major axis, see figure[l] = ujt is 
the orbital phase, and r is the proper time measured by the on-board clock. 
We calculate the gravitational gradients measured along the above orbit for 
the Earth pointing option of S/C attitude. The local Earth pointing frame 
attached to the S/C is determined by the following local tetrad, that we set 
(ei) a along the direction of motion of the S/C, (e 2 ) a along the radial direction 
r, (e 3 ) a = (ei) a x (e 2 )° and (e 0 )° = Z a the 4-velocity of the S/C, see again 
figure{l] Here we do not employ the initial guided option of the S/C attitude 
since the gradient signal from the gravitomagnetic sector will then have twice 
the orbital frequency |29| and can hardly be separated from the far more larger 
signal produced by Earth J 2 component at the same frequency band. Also, as 
one will see, with the help of Earth pointing attitude option, the gradient 
signal of the CS gravity can be separated from that of the semi-conservative 
metric theory. 

The geodesic deviation equation reads 

Z b S7 b Z c \7 c C + R bcd a Z b Z d £ c = 0, (11) 


where £ a is the connection vector determined by the axis of the SGG. The 
gravitational gradient tensor or the tidal tensor is then defined as 

KJ l = R p p X a ZrZ x e,feJ, ( 12 ) 

here we use ~ to mark tensor components under the local tetrad (e^)° to dis¬ 
tinguish the components under the original PN basis (-gf^) a , and e ; /‘ denotes 
the transformation matrix from local frame to the Earth centered PN system, 
and e/' denotes the inverse. From dimensional analysis, we have 


,M, 


, M 2 , 


KJ i = ^[0(—) + 0(v i ) + 0{—) + 0( 
gt a a z 


Mv 2 s , 


Newtonian terms 


Post—Newtonian terms 
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which means that to calculate the 1PN tidal components, that the -4-0(e 4 ) 
terms, only the Keplerian orbits of the S/C with fixed orbital elements are 
needed. To evaluate the Newtonian tidal components, the 1PN correction to 
Keplerian orbits is then needed. However, since the Newtonian components 
depends only on the semi-major a and the unit position vector r, one can then 
write r in terms of the instantaneous orbit elements evaluated at the time of 
observation to work out the Newtonian tidal components. The full solutions 
of the 1PN orbits in CS modified gravity will be the subjects of a separated 
publication. 

The 4-velocity of the S/C reads 

dt 3 3 3 3 

za = Jr ( Q~ t )“ + a “ [- sin ^ Y “ + cos i cos V (^)'“ + sin ■i cos 9 (^)“].. 

(13) 

The ratio can be derived from the line element dr 2 = —g^dx^dx" along 
the circular orbit, to the required order we have 


dt 

dr 


1 + 


0 (e 


(14) 


For the tetrad attached to the reference mass defined in the last subsection, 
we first set (eo) a = Z a , and following the Gram-Schmidt process the three 
spacial tetrad can be worked out to the required order as 


, 9 a 2 ui 2 7 M.. . r , 9 . T . 9 

( e i) = aw (m) +(! + —-—)[- sm!f r (^ r ) +cos JC osf(^) 


'9C 

+ sinicosf(^) a ], 


(e 2 r = 

(1 _ iM. 

V a 

(e 3 ) a = 

/'l _ l' M 


^ o 1 J ^ ^ U 9 , 

3x L 3x z 


, 9 . d , 


" dx 2 ‘ 


y dx 3 ' 


(15) 

“L 

(16) 
(17) 


The transformation matrices e^ a and e,/' are worked out to 1PN level, please 
see Eq.(271 and Eq.(28 ) in Appendix B 

With Eq.( 12), Eq.(131, and Eq.(27 l-(28j, the explicit forms of the 3-dimensional 
1PN tidal matrix in the Earth pointing local frame can be worked out. For 
clarity, we move the explicit expressions in AppendixjB] please see Eq. (29 )- 
(331 for details. These tidal matrices agree with the results in I27IE51E5] that, 
in measuring the off-diagonal components of Kij, the gradient signals from the 
gravitomagnetic sector can be separated from those produced by the Newto¬ 
nian and 1PN gravitoelectric sectors. Now, the most interesting result turns 
out to be that the gradient signature in non-dynamical CS theory can be sepa¬ 
rated from the gravitomagnetic tidal signals of the st and ard parit y-co nserving 
metric theories in K 12 and K 23 components, see Eq. (31) and Eq. (32). 

With a proper orientation of a two-axes diagonal-component SGG, one 
can read out these off-diagonal components produced by the CS modification. 
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Fig. 2 The configuration of the two-axes SGG along the circular orbit. The two perpendic¬ 
ular axes p and q, Eq.( |18[ ), lie within the orbital plane, and they both form a angle of 7r/4 
with the tangent vector of the orbit. 


Denote p and q as the orthogonal normalized 3-vectors of the two axes of the 
on-board SGG, and in the S/C local frame we set 


P W y/? 0} ’ 


q { v / 2’ \/2’ 0} ’ 


(18) 


please see Fig[2] According to the tidal matrices Eq. (29 )-Eq.(33), the sum and 
difference of the SGG read outs I\pp and /\ q q along these two axes can be 
written down as 


K+ = kpp + A qq 


M (4/3 + 2 7 - 3 )M 2 (7 + 2)w 2 


cr 


+ - 


3JcoAcosi 


2 Wn 


3 sin /cos ^ 


K~ = K- - — K~ - = — 

— ^ v pp Jv qq 


0 2 a 3 

Jx'w sin i sin W 
2 a 3 


(19) 

( 20 ) 


here u >0 denote the angular velocity of the rotating local frame relative to the 
parallel transported frame. Please see AppeiidixJU] for the explicit form of the 
readouts along p and q. Therefore, for this simple two-axes SGG settings, the 
gradient signal produced by the CS modification can be isolated, especially 
in K ~, from those produced by the gravitomagnetic sector of the standard 
parity-conserving metric theory. 

For the more reliable three-axes SGG options studied in[27]125][151, we now 
set the three orthogonal SGG axes in the S/C local frame as 


n = 
P = 

q = 


{sin</, —cos (j), 0}, 

/= {cos </, sin (j), -1}, 

/= {cos </, sin (j), 1}, 


( 21 ) 

( 22 ) 


(23) 
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Fig. 3 The configurations of three-axes SGG along the circular orbit. As in Eq.(|2l|-([23j, 
The two perpendicular axes p and q are symmetric to the r — v plane, and n is orthogonal 
to the p — q plane. The angle between n and — r is denoted as 4>. 


where <f> is the angle between n and — r. please see Fig(3| As before, one can 
use the combinations of Kaa, Kf,p and Aqq to single out the interested signals 


K+ = „ a. if ~ . 

— Jv pp ' ^ V qq 

M(3 cos 2<^+l) 2 ((47 — 1) + (8/3 + 87 — 7) cos2<^)M 2 

2 a 3 UJ ° 4 a 4 

((7 + 2) cos2</> + 3 JwA cos i cos 2 </> 

4a a 3 

Joj(x sin i(3 cos 2(t> — 1 ) cos<F — 2 ^ sin i sin 0 cos <p sin A) 

+ " 4a 3 ' ’ 


(24) 


K~ = App - A^q 


3 JujA sin i (cos (j> cos A — 3 sin cf> sin A) Ju>x cos i cos </> 
a 3 + 2 ^ ’ 


(25) 


= A, 


pp + A'qq + /\ n 

J\uj sin i cos A 


M 2 (2/3 — 7 — 1) ( 7 -I )Muj 2 


— 2w 0 . 


(26) 


Please see Appendix|C] for the explicit form of the readouts along n, p and 
q. Therefore, without any change of the optimized three-axes configurations 
studied in [27 | I35 1 I35 ] . the above results turns out to be quite promising. The 
only periodic signal in the difference K~ is from the parity-conserving metric 
theories, which can be distinguished from the DC signal produced by the 
CS modification and then be used to test the gravitomagnetic sector of such 
theories as suggested in |271l25lE51l30] . For the key result here, the only periodic 
signal _ sirw cos & j n sums A'© (or the trace of the total tidal 

matrix ) comes from the parity-violating CS modification, which can be isolated 
from the rest DC terms contained in AT®. 
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3 The measurements of the Chern-Simons gradient signal 


Recovering the SI units, the CS gradient signal to be measured is proportional 


to snyma _ To es ti m ate the magnitudes of the signals and the accuracy 


GJxw sin i cos ^ 

c 2 a 3 

requirements for our experiment, we here assume a polar circular orbit with 
specific altitude as 650km, that a = 7020 km, just like in the GP-B mission. 
For the three-axes SGG, the magnitude of the CS gradient signal is about 
1.5 x 10~ 17 xs _2 with frequency about 0.17 mHz. As reported in (35], the fully 
developed SGG at the University of Maryland with mechanically suspended 
test masses has the sensitivity of 10 -11 s -2 /i7z2 in the signal frequency band. 
For one year lifetime missions, that t ~ ir x 10 7 s, the CS parameter \ can, in 
principle, be constrained to be smaller than 10 2 . From Eq. (j7j , recovering the 
SI units, we have 

»^ 4Ac 

M CS = -• 

Xa 


Therefore, for this conservative option of the SGG sensitivity, the CS mass 
scale can be constrained to be larger than 1.1 x 10~ 15 eU. This is of course a 
rather weak bound on the CS scalar 6 or the mass scale of the theory. While, 
as mentioned before, the high precision SGG with performance level about 
10 ~ li s~ 2 Hz ~2 in the band 10 ~ 4 Hz ~ 10 ~~ 1 Hz has already been developed 
at the University of Maryland E3G2], and even for the performance level of 
10 ~ 15 s” 2 i7z _ 5 at the signal frequency band is within the capability of the 
SGGs under development [35]. To be optimistic, with the likely 4 orders of 
magnitude improvement in the SGG performance level in the future, the CS 
mass scale Mcs could be constrained to be lager than 1.1 x 10 -11 eU for one 
year lifetime missions. 

The signal to be measured will be polluted by various noises, especially 
those from the misalignment between axes and Earth gravity multiples. One 
important point is that the pointing error of the axes will not affect our mea¬ 
surement, since the combination AT® is the trace of the tidal matrix that is 
invariant under axes rotations. For our measurement is similar to that of the 
gravitomagnetic signal in the combination of K~ as discussed in [27 1 135 1 153], 
one can then consult EM for the analysis of the noises from misalignments 
and Earth multiples. As pointed out in m, for Earth pointing orientation, 
the couplings between the alignment error with Newtonian and 1PN gravito- 
electric gradients will not contribute since the alignment does not modulate 
the these gradients in such orientation. While, for the gravitomagnetic gradi¬ 
ent from the parity-conserving theories, its couplings with the alignment error 
will enter into the trace signal A'®. But, to reach the constraint that x < 1CU 2 
mentioned above one only need the misalignments to be measured to lCU'Vad 
in one year, which is rather an trivial task in our experiment. Thus, the test 
of the CS modified gravity will not suffer from the difficulties of meeting the 
stringent alignments and pointing requirements as in the measurements of the 
gravitomagnetic signal of standard parity-conserving theories. The detailed 
analysis concerning all these error sources will be the subjects of a future 
work. 
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At last, in the measurements of gravitomagnetic effects from parity-conserving 
theories, secular signals with increasing magnitude exist [25]|26] due to the 
frame-dragging effect that causes the SGG axes to tilt relative to Earth with 
angular velocity ~ fcmm- These secular signals make the detections of the 
gravitomagnetic effect much easier. The non-dynamical CS modified gravity 
differs with the parity-conserving metric theories only in the gravitomagnetic 
sector. Therefore, it is natural to ask whether or not similar secular gradient 
signals exist for the non-dynamical CS modified gravity, since the axes (or 
gyros) may tilt differently along the orbit due to a different gravitomagnetic 
torque. The answers to this problem require the studies of the 1PN geodesic 
motions and spin precession in the CS modified gravity, which themselves form 
an interesting subject and will left in our future works. 

To conclude here, we studied in this paper the theoretical principles for test¬ 
ing non-dynamical CS modified gravity with SGGs in space. For the two-axes 
and three-axes diagonal-component SGGs, we worked out the characteristic 
signals from the CS modification in the combinations of the measurements of 
different axes, which can be clearly isolated in these combined data. Further¬ 
more, for the three-axes SGG option, the precision tests of the gravitomagnetic 
effect in parity-conserving metric theories[?7][SH][3Hi and the test of the effects 
from CS modification can be carried out independently in the same time. High 
precision multi-axes SGGs are powerful tools for relativistic experiments, es¬ 
pecially that carried in space, and it is promising to add the ability of testing 
the string theory inspired CS modified gravity to the space based experiments 
with orbiting SGGs. 

Acknowledgements This work was supported by the NSFC grands No. 41204051. and 
No. 11305255. 


A The standard PPN metric 

The standard PPN metric has the form m 


goo = -1 + 2 U - 2 PU 2 - 2& w + (2 7 + 2 + o 3 + Cl - 2£)#i 

+2(3y -2/3 + 1 +( 2 + OAl + 2(1 + Cs)^3 + 2(3 7 + 3Cr - 2C)<2> 4 
— (Cl - 2C)_4 — (ai — 02 — a 3 )w 2 U — a 2 w l w :i Uij + (2o 3 - ai)w l Vi + 0(e 6 ), 



- (1 + “2 — Cl + 2 £) Wi 


--(a i — 2a 2 )w l U — a 2 w :i Uij + 0(e 5 ), 
9ij = (1 + 2-yU)&ij + 0(e 4 ), 


1 
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Parameter 

Bound 

Experiment 

7-1 

2.3 x 10 o 

time delay in Cassini tracking 


2 x 10“ 4 

light deflection in VLBI 

p-i 

8 x 10“ 5 

perihelion shift 


2.3 x 10 “ 4 

Nordtvedt effect 

c 

icr^ 

spin precession of millisecond pulsars 

a 1 

4 x icr 5 

orbital polarization of PSR J1738+0333 


10~ 4 

Lunar laser ranging 

0:2 

2 x 10 “ a 

spin precession of millisecond pulsars 

03 

4 x 10“ 2U 

pulsar spin down statistics 

Cl 

0.02 

combined PPN bounds 

C 2 

4 x 10 ~ 5 

binary acceleration of PSR 1913+16 

C3 

Icr 8 

Lunar acceleration 

C4 


not independent 


Table 1 Current values of PPN parameters. 


where the PN potentials read 


f P' j3 / 


r p ' v ' 2 

-J |x-x+ X ’ 

|x- x+ 

f P' U ' d 3 x , 


f P ' n ' d?x‘ 

J |x - x' 

Ix-x'| 

- [, p ' „dV, 

Vi = I 

■, p ' v,i ,dv 


|x-x'| 


|x-x'| 


Wi= , VtvMx-xQ] 
|x — x'l 05 




Uij ~l |x-x '| 3 ** 


p'[v'- (x-x ')] 2 ,3 


= 1 " ^ 


|x-x'|3 |x' - x"| |x - x"| 


)d?x' d?x n 


The matter variables are the rest mass density p , pressure p, coordinate velocity of the 
matter field u*, internal energy per unit mass 77 and the coordinate velocity of the PPN 
coordinate system relative to the mean rest-frame of the universe w z . The PN orders read 


D~c>(e), v 2 ~ U ~ n ~ - ~ C>(e 2 ). 

P 


The standard PN parameters { 7 , (3 , £ , oi, 02 , 03 , Ci> C 2 7 C3 5 C 4 } have the follow¬ 
ing meanings. The parameters 7 and f$ are the usual Eddington-Robertson-Schiff param¬ 
eters used to describe the “classical” tests of GR and are in some sense the most impor¬ 
tant ones. For GR 7 = f3 = 1 are the only non-vanishing parameters. The parameter £ 
measures the preferred-location effects, {ai, 02 , £* 3 } measure the preferred-frame effects 
and{c>: 3 , £ 1 , £ 2 , C3> C 4 } measure the violations of global conservation laws for total mo¬ 
mentum. The up-to-date values of these parameters are summarized in Tab^ |39|. 
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B The transformation matrices and tidal matrices 


Prom Eq. P} -Eq.ljlTJ), the transformation matrices between the local frame and the Earth 
centered PN system up to 1PN level read 


e = 

/! + 


a 2 uj 2 , M 
2 ^ a 


—acj sin d' acj cos z cos d' auj sin z cos d? 

aw — (1 + 0,2 £ 2 — TM) s i n i p (i _(_ a 2 ” 2 _ 7M) cos z cos d/ (1 + a Y — sin z sin 1 


^1 — costf' ^1 — cos z sin ^1 — sin z sin ^ 




_ 7MV 


1 —I cosz 

a I 


e = 


—aoj 0 0 

arnsin!? -(1 + ) sin!? (1 + ^) cosip 0 

—aui cos i cos <P (1 + Sl -! e —f cos j CO s !? (1 + -^) cos i sin >P — (1 + s in j 

y — ausinicos<P (1 + —f ) sin«cos i? (1 + ) sin i sin!? (l + 3^)cosi y 

(28) 
_ r N 

__ _ _ '>)-“ij 

K GE + K GM + K^ s along the circular orbit in the Earth pointing local frame can be 
worked out as 


M 


K N = ^- I 0-2 0 
" 3 ' 0 0 1 


1 0 0 


/ (2/3+37—2 )M 


fP Ad 

k ge = — 

n. 


rr.3 


(6/3+57-5)M 


- (7 + 2)a 2 o ; 2 


(29) 

\ 


( - 2 g-3-y+ 2 )M + (27 + 1)w 2 a 2 y 
(30) 


0 — | A sin z cos ^ \ 

3/1 cosz |zAsinzsin^ 


\ — | A sin z cos d' | Z\ sin z sin d/ —3/1 cosz / 
/ 0 ~\x sin ® sin ^ — IX cos i \ 




Jo; 


— j % sin z sin d/ —^x sin i cos & 


0 


V -ixcosi 


0 


(31) 


(32) 


- Ixsinicosfy 

Here, K GE , K GM and denote the gravitational tidal matrices from the Newto¬ 

nian force, the 1PN gravitoelectric force, the gravitomagnetic force and the contributions 
from the CS modification. Since the Earth pointing frame is rotating relative to parallel 
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transported frames (Fermi shifted), we need to include the tidal matrix from the centrifugal 
force produced by such rotation 


K w 



(33) 


where cjq = + ^0 (e 3 ) denote the angular velocity of the rotation of the local frame. 


C The explicit form of the gradients readouts 

For the two-axes SGG discussed in Subsection |2.3[ see Eq.( |l8[ ) and Fig [ 2 ] the readouts along 
p and q are 


= —0 = - 


M 


PP qq 2a 3 ’ 


f>GE _ £-ge _ (4/3 + 27 - 3 )M 2 _ (7 + 2 )Mui 2 

^pp - qq - 2a 4 2 a 

f>GM _ f>GM _ 3 AJu cos i 
PP qq 2a 3 

f>CS x Jo; sin i(3 cos ^ + sin#') f>cs x^ <jJS ^ rL ^(^ cos ^ r ~ simP r ) 

pp = 4^ 3 ’ ™ = 4 ^ 

For the three-axes SGG in Subsection |2.3| see Eq. ( |21| )-( |23[ ) and Fig. & the readouts along 
n, p and q are 

U ) 2 

K ~ U) 2 T>U) . T>CJ ^0 

nn ^Cb ^PP ^ Y qq 2 ’ 

- N M(3 cos 20 + 1) - N - N M(3cos20 +1) 

K ™ = -2^-’ A PP =A ^= -i^ 3 -’ 

f-GE ((6/3 + 57 — 5) cos 2 0 — (2/3 + 37 — 2 ) sin 2 0) M 2 (7 + 2) cos 2 tf>Mui 2 
= a 4 a ’ 

f>GE = ^.ge = ((47 - 1) + (8/3 + 87 - 7) cos 20)M 2 ((7 + 2) cos 20 + 3^)Muj 2 

pp « 4a 4 4a 

f^GM 3/1 Ja> cos * cos 2 0 


- q m 3AJu (sin j(3 sin 0 sin { P — cos 0cos 'P') + cos / cos 2 0) 

A pp = 2 a 3 ’ 

~ gm 3 AJu: (sini(—Ssin^sin?^ + cos 0 cos + cos i cos 2 </>) 

^ = 2^3 ’ 

~C s J\u> sin i cos 0 (sin 0 sin 'P 3 cos 0 cos S') 

= 2a 3 ’ 

j^CS ( s * n * ((1 — 6 sin 2 0) cos'// — sin 0 cos 0 sin 'P) + cos / cos 0) 


j-^CS J X“ ( s i n / ((1 3 sin 2 0) cos l P — sin 0 cos 0 sin <p) — cos / cos 0) 


4a 3 
— s 
4a 3 


References 

1. M. Niedermaier, M. Reuter, Living Reviews in Relativity 9, 5 (2006). DOI 10.12942/ 
lrr-2006-5 





















Testing CS gravity with orbiting SGG 


15 


2. S. Deser, R. Jackiw, S. Templeton, Annals Phys. 140, 372 (1982) 

3. B.A. Campbell, M.J. Duncan, N. Kaloper, K.A. Olive, Physics Letters B 251, 34 (1990). 
DOI 10.1016/0370- 2693(90)90227-W 

4. B.A. Campbell, M.J. Duncan, N. Kaloper, K.A. Olive, Nuclear Physics B 351, 778 
(1991). DOI 10.1016/S0550-3213(05)80045-8 

5. R. Jackiw, S.Y. Pi, Phys. Rev. D 68, 104012 (2003). DOI 10.1103/PhysRevD.68.104012. 
URL http://link.aps.org/doi/10.1103/PhysRevD.68.104012 

6. S. Alexander, N. Yunes, Phys. Rept. 480, 1 (2009). DOI 10.1016/j.physrep.2009.07.002 

7. T.L. Smith, A.L. Erickcek, R.R. Caldwell, M. Kamionkowski, Phys.Rev. D77, 024015 
(2008). DOI 10.1103/PhysRevD.77.024015 

8. I. Ciufolini, E. Pavlis, Nature 431, 958 (2004). DOI 10.1038/nature03007 

9. I. Ciufolini, E.C. Pavlis, J. Ries, R. Koenig, G. Sindoni, A. Paolozzi, H. Newmayer, 
in Astrophysics and Space Science Library , Astrophysics and Space Science Library , 
vol. 367, ed. by I. Ciufolini, R.A.A. Matzner (2010), Astrophysics and Space Science 
Library , vol. 367, p. 371. DOI 10.1007/978-90-481-3735-0_17 

10. C.W.F. Everitt, D.B. Debra, B.W. Parkinson, J.P. Turneaure, J.W. Conklin, M.I. 
Heifetz, G.M. Keiser, A.S. Silbergleit, T. Holmes, J. Kolodziejczak, M. Al-Meshari, 
J.C. Mester, B. Muhlfelder, V.G. Solomonik, K. Stahl, P.W. Worden, Jr., W. Bencze, 
S. Buchman, B. Clarke, A. Al-Jadaan, H. Al-Jibreen, J. Li, J.A. Lipa, J.M. Lockhart, 
B. Al-Suwaidan, M. Taber, S. Wang, Physical Review Letters 106(22), 221101 (2011). 
DOI 10.1103/PhysRevLett.106.221101 

11. N. Yunes, D.N. Spergel, Phys. Rev. D 80, 042004 (2009). DOI 10.1103/PhysRevD.80. 
042004. URL http://link.aps.org/doi/10.1103/PhysRevD.80.042004 

12. Y. Ali-Haimoud, Phys. Rev. D 83, 124050 (2011). DOI 10.1103/PhysRevD.83.124050. 
URL http://link.aps.org/doi/10.1103/PhysRevD.83.124050 

13. N. Yunes, F. Pretorius, Phys. Rev. D 79, 084043 (2009). DOI 10.1103/PhysRevD.79. 
084043. URL http://link.aps.org/doi/10.1103/PhysRevD.79.084043 

14. Y. Ali-Haimoud, Y. Chen, Phys. Rev. D 84, 124033 (2011). DOI 10.1103/PhysRevD. 
84.124033. URL http: //link. aps. org/doi/10.1103/PhysRevD. 84.124033 

15. K. Yagi, N. Yunes, T. Tanaka, Phys.Rev. D86, 044037 (2012). DOI 10.1103/PhysRevD. 
86.044037 

16. S. Chen, J. Jing, Class.Quant.Grav. 27, 225006 (2010). DOI 10.1088/0264-9381/27/ 
22/225006 

17. K. Yagi, L.C. Stein, N. Yunes, T. Tanaka, Phys. Rev. D 87, 084058 (2013). DOI 10.1103/ 
PhysRevD.87.084058. URL http: //link. aps. org/doi/10.1103/PhysRevD. 87.084058 

18. F. Vincent, Class.Quant.Grav. 31, 025010 (2013). DOI 10.1088/0264-9381/31/2/025010 

19. C.F. Sopuerta, N. Yunes, Phys. Rev. D 80, 064006 (2009). DOI 10.1103/PhysRevD.80. 
064006. URL http://link.aps.org/doi/10.1103/PhysRevD.80.064006 

20. D. Garfinkle, F. Pretorius, N. Yunes, Phys. Rev. D 82, 041501 (2010). DOI 10.1103/ 
PhysRevD.82.041501. URL http: //link. aps. org/doi/10.1103/PhysRevD. 82.041501 

21. P. Pani, V. Cardoso, L. Gualtieri, Phys. Rev. D 83, 104048 (2011). DOI 10.1103/ 
PhysRevD.83.104048. URL http://link.aps.org/doi/10.1103/PhysRevD.83.104048 

22. P. Canizares, J.R. Gair, C.F. Sopuerta, Phys. Rev. D 86, 044010 (2012). DOI 10.1103/ 
PhysRevD.86.044010. URL http: //link. aps. org/doi/10.1103/PhysRevD. 86.044010 

23. V.B. Braginskii, A.G. Polnarev, ZhETF Pis ma Redaktsiiu 31, 444 (1980) 

24. A.G. Polnarev, in Relativity in Celestial Mechanics and Astrometry. High Precision 
Dynamical Theories and Observational Verifications, IAU Symposium , vol. 114, ed. by 
J. Kovalevsky, V.A. Brumberg (1986), IAU Symposium , vol. 114, pp. 401-405 

25. B. Mashhoon, D.S. Theiss, Phys. Rev. Lett. 49, 1542 (1982). DOI 10.1103/PhysRevLett. 
49.1542. URL http://link.aps.org/doi/10.1103/PhysRevLett.49.1542 

26. D.S. Theiss, Physics Letters A 109, 19 (1985). DOI 10.1016/0375-9601(85)90382-2 

27. H.J. Paik, B. Mashhoon, C.M. Will, in Experimental Gravitational Physics , ed. by P.F. 
Michelson, H. En-Ke, G. Pizzella (1988), pp. 229-244 

28. H.J. Paik, Advances in Space Research 9, 41 (1989). DOI 10.1016/0273-1177(89)90006-9 

29. B. Mashhoon, H.J. Paik, C.M. Will, Physical Review D 39, 2825 (1989). DOI 10.1103/ 
PhysRevD.39.2825 

30. H.J. Paik, General Relativity and Gravitation 40, 907 (2008). DOI 10.1007/ 
sl0714-007-0582-4 



16 


L. Qiang and P. Xu 


31. R. Rummel, W. Yi, C. Stummer, Journal of Geodesy 85, 777 (2011). DOI 10.1007/ 
s00190-011-0500-0 

32. M.V. Moody, H.J. Paik, E.R. Canavan, Review of Scientific Instruments 73, 3957 (2002). 
DOI 10.1063/1.1511798 

33. M.V. Moody, H. Paik, K.Y. Venkateswara, P.J. Shirron, M.J. Dipirro, E.R. Canavan, 
S. Han, P. Ditmar, R. Klees, C. Jekeli, C. Shum, AGU Fall Meeting Abstracts p. A786 
( 2010 ) 

34. P.J. Shirron, M.J. Dipirro, E.R. Canavan, H. Paik, M.V. Moody, K.Y. Venkateswara, 
S. Han, P. Ditmar, R. Klees, C. Jekeli, C. Shum, AGU Fall Meeting Abstracts p. A785 
( 2010 ) 

35. H.J. Paik, SQUID Handbook (Wiley, New York, 2006), vol. 2, pp. 545—579 

36. X.Q. Li, M.X. Shao, H. Paik, Y.C. Huang, T.X. Song, X. Bian, General Relativity 
and Gravitation 46(5), 1737 (2014). DOI 10.1007/sl0714-014-1737-8. URL http: 
//dx.doi.org/10.1007/sl0714-014-1737-8 

37. S. Alexander, N. Yunes, Phys. Rev. Lett. 99, 241101 (2007). DOI 10.1103/PhysRevLett. 
99.241101. URL http: //link. aps. org/doi/10.1103/PhysRevLett .99.241101 

38. S. Alexander, N. Yunes, Phys. Rev. D 75, 124022 (2007). DOI 10.1103/PhysRevD.75. 
124022. URL http://link.aps.org/doi/10.1103/PhysRevD.75.124022 

39. C.M. Will, Living Reviews in Relativity 17(4) (2014). DOI 10.12942/lrr-2014-4. URL 
http://www.livingreviews.org/lrr-2014-4 

40. C.M. Will, J. Nordtvedt, Kenneth, Astrophys.J. 177, 757 (1972) 

41. C. Will, Theory and experiment in gravitational physics (Cambridge University Press, 
1993) 

42. L.I. Schiff, Proceedings of the National Academy of Science 46, 871 (1960). DOI 
10.1073/pnas,46.6.871 



